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We present the complete set of analytical solutions of the geodesic equation in Taub-NUT space- 
times in terms of the Weierstrafi elliptic functions. We systematically study the underlying polyno- 
mials and characterize the motion of test particles by its zeros. Since the presence of the "Misner 
string" in the Taub-NUT metric has led to different interpretations, we consider these in terms of the 
geodesies of the space-time. In particular, we address the geodesic incompleteness at the horizons 
discussed by Misner and Taub [j, and the analytic extension of Miller, Kruskal and Godfrey @], 
and compare with the Reissner-Nordstrom space-time. 
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I. INTRODUCTION 

The Taub-NUT solution of the vacuum Einstein equations exhibits intriguing properties. The double name of the 
solution originates from the two parts the solution consists of, and at the same time highlights its strange features. 
In 1951 Taub Q derived a solution valid in the "inner" region (A r < 0, see Eq. (JTJ), which was interpreted as 
a cosmological model (see e.g. Brill Q and Wheeler Q). In 1963 independently Newman, Unti and Tamburino 
(NUT) [| obtained a metric valid in the "outer" region (A r > 0) which reduces to the Schwarzschild metric in the 
limit of vanishing NUT parameter. Misner Q then suggested to consider the Taub space-time joined analytically to 
the NUT space-time as a single Taub-NUT space-time. 

The gravitomagnetic properties of the Taub-NUT space-time are determined by the nondiagonal part of the metric, 
which features the NUT parameter n as a gravitomagnetic charge. This nondiagonal term implicates a singularity 
on the half-axis $ — n, termed "Misner string" which is distinct from the ordinary coordinate singularity associated 
with the use of spherical coordinates. Misner |7J suggested to evade this singularity, by introducing a periodic time 
coordinate and two coordinate patches. The first coordinate patch covers the northern hemisphere and the singularity 
is located along the axis d — ir, while the second patch covers the southern hemisphere with the singularity extending 
along the axis i? = 0. (This construction is analogous to the way in which the "Dirac string" singularity in the vector 
potential of the Abelian magnetic monopole is eliminated.) A periodic time coordinate is thus the price to be paid 
for a Taub-NUT space-time free of these axial singularities. 

Indeed, this periodic identification of the time coordinate makes the solution look rather problematic for physical 
applications because of the resulting causality violations. Noting that every observer at rest in the coordinate system 
would move on a closed time-like world-line, Bonnor |8[ was led to suggest an alternative interpretation of the metric, 
in order to evade this unsatisfactory property of the Misner interpretation of the Taub-NUT space-time. In his 
approach Bonnor retained the singularity at d = n, and endowed it with a physical meaning: he interpreted it as a 
semi-infinite massless rotating rod. In his interpretation the NUT space-time is then considered as being created by 
a spherically symmetric body of mass M at the origin and by a source of pure angular momentum (associated with 
the NUT parameter) which is uniformly distributed along the ■& = 7r-axis. Thus the "Misner string" is representing 
a physical singularity. Moreover, close to the singular semi-axis there looms a region of space-time containing closed 
time-like curves. Consequently, to obtain a physically more reasonable interpretation of the Taub-NUT solution, 
Bonnor was forced to restrict the coordinate range, as to exclude the singular regions along with the region with 
closed time-like curves. Containing holes, the resulting manifold is then incomplete in the sence that not all geodesies 
can be extended to arbitrarily large values of their affine parameters Q . 

Recently, Manko and Ruiz Q reconsidered the interpretation of Bonnor for the general family of Taub-NUT 
solutions, containing a constant C in the nondiagonal part of the metric, gt v oc 2n(cos$ + C). This constant then 
defines the position of the singularity on the axis of the respective Taub-NUT space-time. For the special cases, 
when C = 1 or C = —1 (as considered in [(| Q) there is only a single semi-infinite singularity on the upper or 
lower part of the symmetry axis, respectively. All other C-values correspond to Taub-NUT solutions with two semi- 
infinite singularities. The general source associated with these singularities is then interpreted as two semi-infinite 
counter-rotating rods (or only a single semi-infinite rotating rod in the special cases |C| = 1) with a finite rod (also 
rotating in general) inbetween them. The total mass of the NUT solutions is M and does not depend on the choice 
of C . In particular, none of the rods yields a divergent contribution to the total mass. For the choice C = 0, the 
angular momentum of the middle rod J2 = Cn\jM 2 + n 2 vanishes, and the diverging angular momenta J\ and J2 
of the semi-infinite rods cancel each other. The analysis of Manko and Ruiz is based on the exact calculation of the 
Komar integrals for the mass and the angular momentum of the three rods constituting the NUT source and their 
total values. The unattractive properties of the Bonnor interpretation are thus retained for this general family of 
Taub-NUT solutions: the axial singularity and the closed time-like curves in its vicinity. 

Presuming the periodic identification of the time coordinate, Misner and Taub [l[ explored the geodesies of the Taub- 
NUT space-time and detected the presence of incomplete geodesies. This may seem surprising, since the origin and 
the axis are regular in their interpretation of the space-time. However, employing Eddington-Finkelstein coordinates, 
Misner and Taub found a singular behaviour of some geodesies at the horizons: as one of the horizons is approached 
for the second time, the modified time coordinate diverges and the affine parameter terminates, thus making these 
geodesies incomplete. Now the same type of behaviour is found for geodesies in the Reissner-Nordstrom space-time, 
and is dealt with by analytical continuation of the Reissner-Nordstrom space-time at the horizons, which consequently 
allows the continuation of these respective geodesies. The obvious question following from this observation is therefore 
whether the Taub-NUT space-time can be analytically continued in an analogous fashion as the Reissner-Nordstrom 
space-time. 

Indeed, in 1971 Miller, Kruskal and Godfrey 0] applied Kruskal-like coordinates and provided an analytic extension 
of the Taub-NUT metric for the entire range of r-values. Consequently, in their approach the geodesies, which 
previously terminated at a horizon, could in principle be extended in these new coordinates. They would no longer be 
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incomplete. However, the analytic extension introduced in [2[ came at the prize that the periodic identification of the 
time coordinate was no longer possible, since it would violate the Hausdorff property of the manifold. Consequently, 
without the periodic identification the singularities on the axis are not evaded, and therefore the extension of Miller, 
Kruskal and Godfrey Q is geodesically incomplete, as they note. 

Leaving aside the incompleteness of certain classes of geodesies, which apparently cannot be avoided in either 
interpretation of the vacuum Taub-NUT space-time, let us now address some further interesting aspects of this 
space-time. First of all we note, that the geodesies of the Taub-NUT space-time share many of the properties of 
the trajectories of charged particles in the field of a ma gne tic monopole. A thorough discussion and comparison 
of these orbits can be found in Zimmerman and Shahir [lj|. One particular such prop erty is, that test particles 
always move on a cone. This was also pointed out by Lynden-Bell and Nouri-Zonos 11], who studied a Newtonian 
analogue of monopole space-times and discussed the observational possibilities for (gravito)magnetic monopoles. In 
fact, the spectra of supernovae, quasars or AGN might be good candidates to infer the existence of (gravito)magnetic 
monopoles. 

On the other hand, the lack of observational evidence for the existence of gravitomagnetic masses might be ex- 
plained by its presumably very large value fl2l ]. and by not sufficiently sensitive instruments jl3j |. Nouri-Zonos and 
Lynden-Bell |14j | suggested to look for gravitomagnetic masses by means of gravitational lensing. By considering the 
propagation of light in Taub-NUT space-time they showed that the NUT deflector influences the gravitational lensing 
effect on the null geodesies. In contrast to the lensing in Schwarzschild space-times, the gravitomagnetic field shears 
the observed form of the source. The possibility to detect and the method to measure the gravitomagnetic masses 
with the next generation of microlensing experiments was studied in [l5j . In (l6j even the anomalous acceleration of 
the Pioneer spacecraft was associated with a gravitomagnetic charge. 

We present in this paper the complete set of analytic solutions of the geodesic equation in Taub-NUT space- 
times in terms of the Weierstrafi elliptic functions (the complete set of analytical solutions of the geodesic equation 
in Plebanski-Demianski space-times which Taub-NUT is a special case of, can be found in [l7[). We classify the 
orbits by means of the analysis of the underlying polynomials which include the parameters of the metric and the 
physical quantities characterizing the test particles. In particular, we discuss the incomplete geodesies, arising in 
Eddington-Finkelstein coordinates, as one of the horizons is approached for the second time, and we address the 
analytic extension of the space-time necessary to extend the geodesies, but forcing us to retain the singularity on the 
axis. We also discuss observable effects in Taub-NUT space-times on the basis of analytically defined observables, 
thus making other approaches to experimental signals of the NUT parameter, e.g. [18l Il9|. more rigorous. 



II. THE GEODESIC EQUATION 

The general Taub-NUT solution of the Einstein field equations is described by the metric 

A o 2 
ds 2 = -£(dt- 2n(cosi? + C)dipf - ^—dr 2 - p 2 (dd 2 + sin 2 ddip 2 ) , (1) 

P £\r 

where p 2 = r 2 + n 2 and A r = r 2 — 2Mr — n 2 . Here M is the (gravitoelectric) mass of the solution, and n is the NUT 
charge, regarded as gravitomagnetic mass. For n^O there are always two horizons, defined by A r = 0, and given by 



r± = M ± VM 2 + n 2 . (2) 

One horizon is located at r + > 2M, the other a10 — \n\ < r_ < 0. Between the horizons the radial coordinate r 
becomes timelike, and the time coordinate t spacelike. 

For test particles or light (g^u^u" = 6, 6 = 1 for massive test particles and 6 = for light) moving on a geodesic 
one immediately infers the conserved energy E and angular momentum L 



p 2 \ dr dr 



i A r f dt n ^ r a ( ^dip\ 2 ■ 2 a dip 



L = 2n(cosd + C)-± — - 2n(costf + C)-f- +p 2 sin 2 ^-^ , (4) 
p z \dr dr J dr 



where r is an affine parameter along the geodesic. 



1 Prom r_ = M — \J M 2 + n 2 we obtain r 2 _ — 2r_M = n 2 . Since r_ < it follows that — \n\ < r_. 
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For convenience, we introduce dimensionless quantities (rs = 2M) 

r = — t = — f = — n = — L = — (5) 

rs r s rs ' r s rs 

It is straightforward to see that in Taub-NUT space-times the Hamilton-Jacobi equation 

is separable and yields for each coordinate a corresponding differential equation 

= R (7) 



2 



= e (8) 

dip V — 2hE cos i? 



d"f sin 2 ■& 



, with L' = L — 2nEC (9) 



— = EKr- +2n(cosi? + C) ^ . (10 

d 7 A r sin 2 z9 v ; 

Here we used p 2 = r 2 + h 2 and A r = r 2 — r — n 2 and introduced the Mino time 7 through p 2 dj = dr [20]. We also 
defined 

R= (r 2 + h 2 ) 2 E 2 - A r (6 r 2 + L 2 + k) (11) 

e = k-Sn 2 + L 2 -( L/ - 2 ^ COS ^ . (12) 

sin 2 V ; 

The separation constant k is known as Carter constant. A study of the geodesic equation in Weyl coordinates has 
been performed in [21] . 

A. Gauge transformation 

Consider a transformation t = t' + 2nCip. It brings the metric ([T]) into C-independent form: 

A n 2 
ds 2 = -£ (dt' - 2ncas<&dtp) - ^—dr 2 - p 2 (dd 2 + sin 2 <&d<p 2 ) , (13) 

P L±r 

with corresponding conserved energy E' = E and angular momentum LI = L — 2hEC for a test particle or light. Thus, 
the apt notation L' introduced in the Hamilton- Jacobi equations (l7])- (jTU]) represents a conserved angular momentum 
in the gauge transformed metric (|13[) . 

The Hamilton- Jacobi equation for each coordinate of the metric (fT3"|) yields 

dr\ 2 

f) =e 



) =R (14) 

2 

dj / 

dip L'-2nEcos'd ,„„. 

— = t> (16) 

fry sin v 

dt' p 4 V - 2nE cos -d , x 

— = E^r- + 2n cos 5 , 17) 

d-y A,. sin 2 ?? V ; 

with 

R = {r 2 + n 2 ) 2 E 2 - h r {5? 2 + L' 2 + k') (18) 

G^k'-Sn 2 + L' 2 -^-^p^l. (19) 

sin v 
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The Carter constant k' is related to k as k' — k = L 2 — L' 2 . 

From the considerations above one concludes that the metric ((TJ can be brought into form (|T5)) which corresponds 
to the metric ((TJ) for C = 0. The same feature is observed for the Hamilton- Jacobi equations (fT))- (fTU)) and (fTl)) - (fTT)) 
under interchange of the constants L' and L, k' and k and the coordinates t' and t. 



III. COMPLETE CLASSIFICATION OF GEODESICS 



Consider the Hamilton- Jacobi equations ([7))- (ITUl) . The properties of the orbits are given by the polynomial R (fTTj) 
and the function O (|12[) . The constants of motion (energy, angular momentum and separation constant) as well as 
the parameters of the metric (dimensionless NUT-parameter) characterize these polynomials and, as a consequence, 
the types of orbits. In this section we discuss the motion in Taub-NUT space-times in terms of the properties of the 
underlying polynomial R and the function 0. 



The ■&— motion 



In order to obtain from Eq. (JSJ) real values of the coordinate •§ we have to require O > 0. This implies c 2 := 
k — Sh 2 + L 2 > 0. With the new variable £ := cosi9, Eq. ([5]) turns into the equation 



rff 2 

ch, 



5 with C := < 2 + b£ + c , (20) 



with a simple polynomial of second order on the right hand side, where a = — (c 2 + Ah 2 E 2 ), b = AnEL', and 
c — c 2 — L' 2 . Since c 2 > we have a < 0. This means that O^ can be positive if and only if there are real zeros of 0£. 
The polynomial 0^ plays the role of an effective potential for the ^-motion. The zeros of define the angles of two 
cones which confine the motion of the test particles. (A similar feature appears in Kerr space-times [HI].) Moreover, 
every trajectory is not only constrained by these cones but lies itself on a cone in 3-space [H, H3, El (the discussion 
for C = —1 in the papers [l0( can be extended for the general family of Taub-NUT solutions). 

The discriminant D = b 2 — 4ac of the polynomial 0^ can be written as D — 4ciC2 with c\ := c-2 — L' 2 + 4h 2 E 2 . 
The existence of real zeros of 0£ requires D > 0. This implies that both c\ and C2 should be non-negative 

ci = c 2 - L' 2 + Afi 2 E 2 >0 

c 2 = k - 5n 2 + L 2 > . [ ' 

These are conditions on the parameters E, L, and k for some given n. As long as k — Sn 2 is positive there are no 
constraints on L and E. When k — Sn 2 becomes negative the inequalities (|2ip imply lower limits for the energy and 

angular momentum given by L m j n = ±\/8n 2 — k and E m i n = 2 n(i-c' i ) (^~LC ± %/ Sn 2 — k + C2C 2 ^j . 

One can show that 

(L'cosi9-2Eii) 2 (L 1 - 2Encosti) 2 

© = ci r^~ q = c 2 r^~ q ■ (22) 

sin i) sin v 

The zeros of 0£ are given by 



2EhL' T , . 

51.2 - 1 I r,2 ■ \ 26 ) 



Cl 



The function 0^ describes a parabola with the maximum at ( j — pj^, ^ ) ■ For nonvanishing h, E, and L 



The conditions ([2T]) ensure the compatibility of £ 6 [— 1, 1] with 0^ > 0. 

2hEL' c\C2 
^a + L' 2 ' Cl + U 

the maximum of the parabola is no longer located at £ = or, equivalently, the zeros are no longer symmetric with 
respect to £ = 0. Only for vanishing h, E, or L both cones are symmetric with respect to the equatorial plane. 
The i9-motion can be classified according to the sign of c 2 — L' 2 : 

1. If C2 — L' 2 < then 0^ has 2 positive zeros for L'En > and d g (0,7r/2), so that the particle moves above the 
equatorial plane without crossing it. If L 1 En < then $ € (k/2, tt). 



6 



2. If c 2 - L' 2 = then 6 5 has two zeros: = and £ 2 = tt^tE^W- If L '^" ^ then £ G I ' X ) and ® e (°> f ]• 
If L'En < then £ € (—1, 0] and $ € [f , 7r). If L' = 2Eh then the ^-motion fills the whole upper hemisphere 
d G [0, ^]. The motion fills the whole lower hemisphere with d £ [? ,tt] if L' = —2Eh. 

3. If c 2 — L' 2 > then O^ has a positive and a negative zero and $ G (0, 7r), and the particle crosses the equatorial 
plane during its motion. 

In general, the second term of the function O in (|22|) diverges for d — >• 0, 7T. However, if £' = 2£'n this term is 
regular for ■& — and if L' = —2En it is regular for d — ir. J£ L' = ±2En, then ci = c 2 . The regularity of O in these 
cases can be seen from 

e = c 2 -4E^^^l, (24) 
sin v 

by application of L'Hopital's rule. If, furthermore, E 2 = then = k — 5n 2 + L 2 which is independent of *&. 

In the special cases when one of the constants c\ or c 2 or both vanish, 0^ has a double root which is the only 
possible value for i?. One can distinguish three cases: 

1. If d = and c 2 > 0, then £ = for i' 2 > 4E 2 ?i 2 . 

2. If c 2 = and c x > then £ = ^ for L' 2 < AE 2 n 2 . 

3. If ci = c 2 = then £ = ±1 implying that § = or i? = 7r are possible. In this case L' = ±2Eh (as discussed 
above). 

This means that during a test particle's motion the coordinate $ is constant and the trajectory lies on a cone around 
the i? = 0, 7r-axis with the opening angle arccos £. In this case we immediately obtain from (j9|) that ip(j) = C (7 — 7i n ) 
with a constant C = L 7 2 ?F^ which in case 1 is C = V and in case 2 is C = 0. 

Thus, the nonvanishing of ci and c 2 indicates that the motion of the particle is not symmetric with respect to the 
{) = 0, 7r-axis. Therefore, these two constants may be regarded to play the role of a generalized Carter constant which 
appears, e.g., in the motion of particles in a Kerr space-time. 

B. The r— motion 

1. Possible types of orbits 

Before discussing the f-motion we introduce a list of all possible orbits: 

1. Transit orbits (TO) with r G (±00, ^foo) with a single transit of r = 0. 

2. Escape orbits (EO) with ranges 

(a) (r 1; oo) with r\ > or 

(b) (-00, ri) with r[ < 0. 

These escape orbits do not cross r = 0. 

3. Crossover escape orbits (CEO) with ranges 

(a) (ri,oo) with r\ < or 

(b) (-00, ri) with r[ > 0, 

which cross r — two times. 

4. Bound orbits (BO) with range f £ (7*1, r 2 ) with r\ < r 2 and 

(a) either n , r 2 > or 

(b) ri,r 2 < 0. 

5. Crossover bound orbits (CBO) with range r £ (ri,r 2 ) where n < and r 2 > 0. 

The mathematical condition from the geodesic equation for crossing the origin r = is R(0) > 0. This implies 
i?(0) = h 2 (n 2 E 2 + c 2 + Sh 2 ) > 0, which is always fulfilled under the conditions (f2~Tj). 
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2. The radial motion 



The right hand side of the differential equation ([7]) has the form R = 5Z i=0 bif 1 with the coefficients 

b A = E 2 -8 (25) 

b 3 = 8 (26) 

b 2 = 2h 2 E 2 - c 2 (27) 

bi = c 2 + 8fi 2 (28) 

6 = n 2 {n 2 E 2 + c 2 + <5n 2 ) . (29) 

Let us now consider massive particles only, that is 8 = 1. In order to obtain real values for f from ([7]) we have to 
require R > 0. The regions for which R > are bounded by the zeros of i?. The number of zeros depend on the 
values of E, L, and n. For E 2 -1^0 there are regions with j and j ' ± 2 real zeros whose boundaries are given by 
R = and i?' = 0. This is used for the parameter plots shown in Fig. [1] One has to additionally take care of the 
change of the sign of E 2 — 1 when E 2 crosses E 2 — 1. Then the sign of R(f) for r — > ±oo changes. Furthermore, the 
case E 2 = 1 requires additional attention: If the line E 2 — 1 is contained in a region with 4 or 2 zeros then on this 
portion of the line we have one zero less, that is only 3 or 1 zeros, respectively. Taking all these features into account 
we obtain the L-E diagrams of Fig. Q~] 

Before classifying all orbits we show that it is not possible to have an orbit completely constrained to the region 
between the two horizons. If such an orbit would exist then the turning points should lie in this region. These turning 
points are given by 

df\ 2 , o 9\ 2 / 9 x r 2 + L 2 + k\ , „, 

— = if 2 + n 2 ) E 2 - A r »- , 30 

d-yj V ' \ (f 2 + n 2 ) 2 J 

where L 2 + k > as can be inferred from ([2~Tj) . Then it is clear that for E 2 = we have A,. = which gives the 
horizons r = r± as turning points. This case requires k > n 2 . If E 2 > then also A r > 0, implying that the turning 
points arc in the outer regions, that is in r > r + and r < r_. Therefore it is not possible to have geodesic motion 
completely constrained to the region between the two horizons. Moreover, having crossed one of the horizons once, a 
particle is forced to pass the origin r = and then cross the other horizon as well. 

Furthermore, it is not possible to have a turning point at one of the horizons only. If, e.g., r = r_ is a turning 
point, then A r (r_) = 0. Then we necessarily have E 2 = 0. In this case also r = r+ is a turning point. Therefore it is 
not possible to have a turning point at one of the horizons only. If E 2 = then the turning points at the horizon are 
the only turning points. That is, if one of the turning points is on one of the horizons then a second one is lying on 
the other horizon, and there are no further turning points. 

The expression 

~ f 2 + L 2 + k . . 

V cS :=A r 2 (31) 

(r z + n z ) 

may be regarded as an effective potential. Though it does not determine the motion of the particles in the usual sense 
it determines the turning points through the condition E 2 — V c r — 0. 

Using the L-E diagrams in Fig. Q] as well as the above considerations we can give all possible combinations of zeros 
of R and an interpretation in terms of specific types of orbits which are summarized in Table [I] Note that at the 
moment we are not discussing features related to singularities and/or geodesic incompleteness. This discussion will 
follow in Section [V] In the present section we are just exploring the types of orbits that are mathematically possible 
as solutions of the geodesic equations. 

The types of orbits related to various parameters are given by: 

• Region (0): no real zeros. We have TOs only with particles moving from ±oo to =Foo which we call orbit A. 

• Region (1): one real zero. Here the coefficient of the highest power in R, given by E 2 — 1 for massive test 
particles, vanishes for E 2 = 1. On the E 2 = 1 line within the region (2) there is one real zero. The orbit is a 
CEO with particles coming from r = +oo. 

• Region (2): two real zeros. Due to the fact that the coefficient of the f 4 -term is given by E 2 — 1 the sign and, 
thus, the type of orbit changes for E 2 larger or smaller than 1. For E 2 < 1 only CBOs are possible. 

- Region (2)+: We obtain EOs and CEOs. 
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type 


rsg-ion 








A 

i\ 


l(Y\ 
K u ) 


n 
u 


U II 1 




B 


(!) 





1 1 II ! 


CEO 


C 


(2) + 





2 • • 1 


EO, CEO 


D 


(2)- 


1 




• CBO 


Dn 






H ■ 


f ' CBOo 


E 


(3) 












2 


i • 






• • ' CBO, EO 


F 


(4) + 


2 


2 * * 1 


• • ' EU, CBO, EO 


G 


(4)- 


3 


1 • 1 


• • • ' CBO, BO 



TABLE I: Types of polynomials and orbits in the Taub— NUT space-time. The thick lines represent the range of orbits. Turning points arc shown 
by thick dots. The regions (1) and (3) are the line E 2 — 1 lying in the regions (2) and (4). The indices + and — indicate E 2 > 1 and 
E < 1, respectively. The position r — is shown by a vertical line, the horizons arc indicated by a vertical double line. In the case Do 

the turning points lie on the horizons for E 2 — 0. 



— Region (2)_: Due to the sign change mentioned above the types of orbits change leading to one CBO. For 
< E 2 < 1 the two turning points are larger than r + and smaller than r_, respectively. For E 2 = the 
two turning points are lying on the two horizons. 

• Region (3): three real zeros. Again, for E 2 = 1 the term of highest power in R vanishes giving three zeros 
for the E 2 = 1-line contained within the region (4). For E 2 > the turning points cannot coincide with the 
horizons. 

• Region (4): four real zeros. Here again the types of orbits change depending on the sign of the coefficient E 2 — 1 
of the leading term in R. 

— Region (4) + : We obtain two EOs and a CBO. For the CBO the turning points cannot be on the horizons. 

— Region (4)_: Here we find BOs and CBOs. The turning points cannot be on the horizons. 

The corresponding effective potentials are exhibited in Figure [5] We note the asymmetry of the possible types of 
orbits with respect to the outer regions. In the negative r region no bound orbits are possible, since the effective 
potential is repulsive here. 



IV. SOLUTION OF THE GEODESIC EQUATION 

Now we present the analytical solutions of the differential equations (|7|)- p0l) . 



A. Solution of the ■&— equation 

The solution of Eq. (|20p with a < and D > is given by the elementary function 

0( 7 ) = arccos(i- (V^sin (v^7 - 7&) " &)) , (32) 



where 7 




and 7; n is the initial value of 7. 



B. Solution of the f— equation 

For timelike geodesies the polynomial R in (fTTj) is of fourth order. A standard substitution r = ±j + ?r, where fa 
is a zero of R, reduces (O to a differential equation with a third order polynomial = R3, where i? 3 = J2i=o biX 1 . 
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(i) n = 0.5, fc = 15.0 



FIG. 1: Parametric L-_E digrams showing the location of the regions (0)— (4), which reflect the number of zeros of the polynomial R in Eq. jilt . 

Each region contains a set of orbits pe culiar to it which are describe d in Tableland in the text (Section |III B 2) . Here the influence of 
the independent var iation of n (fig |l(a){f{d^ , and k (fig |l(c)|[(T)| on the zeros of R is presented. The dashed region denotes the 
prohibited values of E and L following from the inequalities ]2H . Here C — 0. 



A further substitution x = (4y — -y) transforms that into the standard WeierstraB form 



fj) =V-5a»-fl3:=i%(!/), (33) 



where 



_ &2 &1&3 &0&3 6| ,„,v 

■ 92 -12~~' ff3 = ^^~Tcf~216 • (34) 

The differential equation (|33|) is of elliptic type and is solved by the Weierstrafi p-function [23[ 

y{i) = p (7- iLi 52,53) , (35) 

where 7,' = 7 in + f°° , dy with y; n = ±% (f; n — r^) - + f§. Then the solution of J7]) acquires the form 

f = ± JT+rR- (36) 

4p (7-7i'„;52,53) ~ f 
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(a) 



CO 



D 



V otf for h = 0.5, k = 1.0, L = 3.0. 



V off for n = 0.5, k = 1.0, L = 2.0. 





A 




C 



















FIG. 2: V off for orbits of type A, B, C, D, D , E, F, G from Tabic [I] showing the position of the orbit types CBO, CEO, BO, EO, TO. The points 
indicate the turning points of the motion and the dashed lines correspond to the values of the energy E 2 . At infinity the effective 

potential tends to the value lim^^^oo V e ff — 1. 



C. Solution of the 93— equation 

Eq. © can be simplified by using ([5]) and by performing the substitution £ = cos d 

dip = 



where 0£ is given in ([20"]) . This equation can be easily integrated and the solution for a < and D > is given by 



(37) 



I--L 



e(7) 

tin 



(38) 
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where 



It can be shown that 



V ± 2Eh 

1+ = — — - arctan ■ 

± \U ± 2En\ 

2aS, + b 



1-uB 



\ZT=u*Jb% - 1 



B± = 



'D 
b±2a 



D 



B±-l = 



4(-g) 

D 



(2Eh =F L') 2 > 0. 



For the special case c 2 = L and L' = 2Eh the solution reduces to the simple form 

«7) 



y(7 ) = 1 arctan ^ ~ ^ 



(39) 

(40) 
(41) 

(42) 
(43) 



Eq. (ITOl) consists of f and $ parts: 

t — t\ n ~ E I — = 

A r VR 



D. Solution of the t— equation 



, . L'- 2fiEcosd dd 

2n \ (costf + C) 



sin z i? ve 

Again with the substitution £ = cos$ the integral /^(t?) can be easily solved: 



( 7 ) = fi(l+ + /_) |* (7) +Cn(/_ - /+) 



«7) 4n 2 £ . 2a£ + 6 
H ■ arcsm — -== — 



f(7) 

Sin 



with 7± as in ((39 



We now consider If. The substitution f = ±- 



"4y- 



+ ffl reexpresses If in terms of y: 



(44) 



(45) 



7,(7) = E f V T- 



dy 



(f R (4y-%)±6 3 ) +n 2 (4y-f)' 



(46) 



where A y = A r (f fl .) (4» - b ff ± (2f fl - 1) 6 3 (4y - f ) + b\ = 16A r (f R )(y - Pl )(y ~ p 2 ). Here, Pl and p 2 are two 
zeros of A y . 

We next apply a partial fractions decomposition upon Eq. (|46"1) 



If{l)=El T-_ 7 



=J»-Pj fa-Pa)* 



(47) 



where P3 = f§ and -K*, « = 0, ..,4, are constants which arise from the partial fractions decomposition. These 
depend on the parameters of the metric and the test particle and on Pr.. After the substitution y — p(v) with 
p'(v) — (— l) p -\/4p 3 (w) — g 2 p{v) — 93, where p is either or 1 depending on the sign of p'{v) in the considered 
interval and on the branch of the square root, Eq. (|47|) simplifies to 



J vi« \ 4=1 



K; 



K 4 



^p{v)-pj (p(v)~ P3 ) 2 



dv 



(48) 



Here v = v(<y) =7- 7i ' n and v m = v(y m ). 
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(a) 3d plot. The sphere shows the horizon r+. 



(b) Projection onto the x — z— plane. The two 
circles denote r+ and r_ . 



FIG. 3: TO with parameters n = 1, k = 1, L = 2, E — 3. The figures combine both regions r > (blue) and r < (sea green). 



The final solution takes the form (details can be found in the appendices [B] and [Cj) 



Hi) = T(-l) p Et (Ko + A 2 K±) (v - v in ) + 

{ i=i 



£;7 



3 = 1 



P'(Vji) 



C(Vji)(v - V in ) + log 



cr(v m - v^) 



— -2 C(f - ««) - C(«in - «S<) + — r log — 

(p'(V3i)) V p {V3i) CT Win " «3i) 



(49) 



where are the poles of the functions (p(u) — Pj) 1 and (p(u) — in (|47|) such that ^(uji) = = p(^j2) 

since p(v) is an even elliptic function of order two which assumes every value in the fundamental parallelogram 
with multiplicity two. Here £(u) is the Weierstrafi zeta function and a(v) is the Weierstrafi sigma function [23| ; 

P(v3i) . p" \v 3i )((v 3t ) 



A 2 = - V — 

^ \ (trJ 



( v 3i)f (p'(V3i)) 3 



is a constant. 



E. The complete orbits 



With these analytical results we have found the full set of orbits for massive point particles: A TO is shown in 
Fig. El a BO in Fig. g| a CBO in Fig. an EO in Fig. and a CEO i n Fig. The degenerate case of a CBO 
where the turning points are lying on the horizons is exhibited in Fig. 8(a) A further degenerate case, an orbit which 
crosses the $ = 0-axis, is shown in Fig. |8(b)] The CBO, CEO and TO contain a positive r and a negative r par t. For 
simplicity and since it does not create confusion, these two parts are included in the same plot. In Fig. [31 Fig. 
Fig. [7(b)] and Fig. ~ 



5(b)] 



ways 



a) the part of the orbit with negative r is plotted in sea green, whereas the blue color a ~ 

represents positive r. 

The constant C appearing in the general family of Taub-NUT solutions (eq. (p}) is chosen zero in the calculation of 
the geodesies. In general, it follows from the Hamilton- Jacobi equations (JZJ)-® that the constant C can be absorbed 
in the constant L' and does not influence the characteristic features of the geodesies in the Taub-NUT space-times 
discussed in section lHll Moreover, as it is shown in the section lll Al one can switch between the C-dependent metric ([TJ) 
and Hamilton- Jacobi equations (fTf- fTU)) and the C-independent metric (fT5)) and Hamilton- Jacobi equations (|14l) - (fT7|) 
by introducing a gauge transformed coordinate if and interchanging the constants of motion L' and L, k' and k. 

These orbits are just solutions of the geodesic equation for the spatial coordinates. The complete discussion of these 
orbits is possible only by including the solutions for the time coordinate and the proper time, which is critical for the 
issue of geodesic incompleteness and other singularities. This is done in the next section. 
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(a) 3d plot. The sphere shows the horizon r_|_. (b) Projection onto the x — z— plane. The two 

circles denote r_j- and r_ . 



FIG. 4: BO with parameters n = 0.5, k = 1.0, L = 2.0, E 2 = 0.94084. The sphere denotes the horizon at r + . 



3-1 




(a) 3d plot. The sphere shows the horizon r_|_. (b) Projection onto the x — z— plane. The two 

circles denote r+ and r_ . 



FIG. 5: CBO with parameters n — 0.5, k — 1.0, L — 2.0, E 2 — 0.94084. In the 3d— plot the orbit for positive r (blue) is shown. The transition to 

negative r (sea green) can be seen in the projection onto the x — 2-planc. 

V. THE ISSUE OF SINGULARITIES 

In this section we discuss the nature of the singularities which are present in the Taub-NUT space-times. We 
address these singularities in terms of the analytic solutions of the geodesic equation. Different interpretations of the 
Taub-NUT space-time give rise to different types of singularities. 

1. The family of Taub-NUT space-times as given by the metric (QJ possesses a singularity along the whole $ = 0, n 
axis when |C| 7^ 1 [9j, which includes our choice C = 0, Eq. (Jl]). The choice |C| = 1 leads only to singular 
semi-axes (see e.g. Misner 0]). 

2. In the interpretation of Misner Q and Misner and Taub [l[ one can avoid the singularity present on the axis by 
employing a periodic identification of the time coordinate. Hypersurfaces r = const then possess the topology 
S 3 . In this case it is, however, not possible to eliminate the geodesic incompleteness which arises at the horizons, 
as discussed below. 

3. When analytically extending the space-time at the horizons by gluing further copies of the space-time following 
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(a) EO with parameters h — 0.5, k — 1.0, (b) Projection onto the x — z— plane. The two 

L — 3.0, E 2 — 1.369. circles denote r+ and r_. 



FIG. 6: An example for an EO. 




(a) CEO with parameters n — 0.5, k — —0.5, (b) Projection onto the x — z— plane. The 

L = 2.0, E = 1.0000001. two circles denote r + and r_. 



FIG. 7: An example for a CEO. The figures combine both regions r > (blue) and r < (sea green). 

Miller, Kruskal and Godfrey @], one has to retain the singularity along the $ = 0, it axis. Therefore one finds 
the space-times considered by Bonnor Q and Manko and Ruiz Q , where the axial singularity is interpreted as 
a physical singularity. Such a singularity might possibly be avoided, however, if an appropriate interior solution 
were included. 

Before discussing the singularities of Taub-NUT space-times in terms of the behavior of the geodesies, we visualize 
Taub-NUT space-times using Carter-Penrose diagrams. 

A. Singularities in a single copy of Taub-NUT space— time 

The Carter-Penrose diagram for a single copy of Taub-NUT space-time is shown in Fig. |H1 It consists of three 
regions: region I with oo > r > r + > 0, region II with r + > r > r_ and region III with 0>r_>r>— oo. This 
diagram is helpful in discussing the singularity issues. One can consider two versions of it because the geodesic motion 
in regions I and III is different. (Recall that the potential in region III is repulsive.) Thus the motion of a particle 
depends on whether it starts in region I or region III. 
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(a) Degenerate CBO with parameters n — 0.5, 
k = 1.0, L = 3.0, E = 0. The turning points 
coincide with the horizons, where r_ < is 
the inner circle. The motion takes place in a 
plane. The figure combines both regions r > 
(blue) and r < (sea green). 




(b)CBO: n = 0.5, k = h 2 , L = 0.9, E = L/(2n). The red 
point indicates the i9 — 0-axis. In the Misncr 
interpretation the axis is regular and a test particle 
can move across $ — 0; alternatively, when the axis is 
interpreted to represent a physical singularity, the 
orbit must end when it encounters this singularity. 
Details are given in SectionM 



FIG. 8: Two degenerate orbits. A CBO with the horizons as turning points and a CBO crossing the $ — 0, it axis. 



The possible orbits obtained above are depicted in Fig. [TO] From this diagram it is immediately clear that the 
CEOs shown in red, as well as the CBOs shown in blue (and the D n shown in green which represent a special case of 
the CBOs) are geodesically incomplete in a single copy of Taub-NUT space. The CEOs, for instance, start at spatial 
infinity, cross the horizons r + and r_, and are forced to end when they reach the horizon r_ a second time. Therefore, 
as long as there is only a single copy of Taub-NUT space, there is no way for such a geodesic to cross the horizon r_ 
a second time, implying geodesic incompleteness of the space-time. This illustrates the problem in the interpretation 
of Misner and Taub, where the periodic identification of the time coordinate precludes the analytic extension of the 
space-time with further copies. 

B. Singularities in the Misner and Taub version of the Taub-NUT space— time 

The geodesic incompleteness in this version of Taub-NUT space-time has been investigated by Misner and Taub 
[l| . We now explicitly demonstrate this geodesic incompleteness by means of our analytic solution: By introducing 

Eddington-Finkelstein like coordinates ip± through 2nip± — t ± J r -fc-dr, one eliminates the singular behavior of the 
metric ([1]) at the horizons A r = and one obtains two metrics 



ds 2 ^An 2 ^-(dijj- (costf 



- C)d(p) 2 T 4n {dtp - (cos + C]dtp) dr - p 2 (dtf 2 + sin 2 Mip 2 



(50) 



With (7) and ^(7) as in (fH)) the equation for the coordinate ip takes the form 



2n(V>± - Vin) = (t - t in ) ± / j-df 



Ir(l)+Uh)± (f+~log|A r | + i V / l + 4n 2 log 



2f-l 
V4n 2 + 1 



(2f-l) 
V4n 2 + 1 



f(7) 



(51) 



Figure [TT] shows the V'(7) _ dependence which reveals the incompleteness of crossover bound geodesies. Consider 



e.g. '0+ in Fig. 11(a) A test particle starting its motion somewhere in the region A r > can cross the horizons r+ 
and r_ only once; when it approaches r_ for the second time the function ^(7) diverges and prevents the particle from 
leaving the region A r < 0. This behaviour represents the geodesic incompleteness observed by Misner and Taub [l|. 
Also if a particle starts from f = the geodesic is incomplete, since once a particle leaves the region A r < after 
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(a) Taub-NUT space-time I (b) Taub-NUT space-time II 



FIG. 9: Two causal versions Taub-NUT space— time in a Carter— Penrose diagram. Here region I is given by r > r+, region II by T- < r < r+ and 

region III by r < r_. For the Kruskal— like analytic continuation sec Fig. 1 12 1 




(a) Taub-NUT space-time I. The CEO (red) terminates at the (b) Taub-NUT space-time II. The CEO (red) starts at the 
horizon r_ , the CBO (blue) starts at the horizon r+ and horizon r„, the CBO (blue) starts at the horizon r„ and 

terminates at the horizon r_, the orbit Do (green) is a special terminates at the horizon r+. 

case of the CBOs. 



FIG. 10: Topology of orbits in Carter— Penrose diagrams of Taub-NUT space-time. The orbits drawn in black arc standard orbits with infinite 

proper time, the orbits in red, blue, and green arc gcodesically incomplete. 



crossing one of the horizons r_ or r+, it cannot cross this horizon a second time, because ^(7) diverges there. At 
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(a) geodesic incompleteness of Taub-NUT space— time for 
V>+(7) 



(b) geodesic incompleteness of Taub— NUT space— time for 

i/>_(7) 



FIG. 11: Taub-NUT space-time: and »/>_ for n = 0.5, k = 1, L = 2 and E 2 = 0.93. The orbit is a CBO of the type depicted in Fig.|5(a) 



the same time the affine parameter 7 cannot be continued further, which indicates the geodesic incompleteness of the 
Taub-NUT space-time. 



Figure 11(b) for ip- demonstrates the incompleteness of the geodesies at the first pair of horizons. 



Thus, the Eddington-Finkelstein-like transformations eliminate the singular behaviour of the original Schwarzschild- 
like coordinates only at the first crossing of the horizons, but lead to incomplete geodesies at the attempted second 
crossing. 

Because an analytic extension of the space-time with a second copy is not possible without destroying the Hausdorff 
property of the underlying topological space Q , the Taub-NUT space-time with periodically identified time coordinate 
is geodesically incomplete. 



C. Singularities in the Kruskal— like analytic extension of the Taub-NUT space— time 



In order to extend these incomplete geodesies, Miller, Kruskal and Godfrey presented a Kruskal-like analytic 
extension of the Taub-NUT space-time, where an infinite sequence of further copies of Taub-NUT space-time is added. 
This procedure is completely analogous to the case of the Reissner-Nordstrom space-time discussed in Appendix fX] 
The Carter-Penrose diagram of this extension of Taub-NUT space-time is presented in Fig. [TJJ 

In this extended space-time it is possible to continue the CEOs and CBOs. These geodesies are therefore no longer 
forced to terminate at a horizon, when trying to pass it for a second time. This is illustrated in Fig. [13l where these 
orbits are displayed in the Carter-Penrose diagram. 

The price to pay for the continuation of the geodesies at the horizons is, however, that the singular $ = 0,7r-axis 
cannot be eliminated by a periodic identification of the time coordinate. Thus, the axis represents a singularity, where 
geodesies must terminate. In the classification of Ellis and Schmidt l28l l29| th is type of singularity has been termed 
a quasi-regular singularity, since the curvature remains finite (see also [30j, |3l|). Nevertheless, a geodesic ends when 
reaching the singularity, analogous to the more often encountered conical singularity. A further discussion of this type 
of singularity may be found in [32[ and [H[ , classifying this singularity further as a torsion singularity. 

Thus, the introduction of the Kruskal-like analytical extension of the space-time eliminates the geodesic incom- 
pleteness at the horizons but only trades it for a physical singularity along the <& = 0, 7r-axis, where the geodesies 
are forced to terminate: consequently the geodesic incompleteness of Taub-NUT space-times is retained in cither 
interpretation of this intriguing vacuum solution. 
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FIG. 12: Kruskal— like analytic continuation of Taub— NUT space— time in a Carter— Penrose diagram. Here region I is given by r > r-f, region II by 
T- < r < r_|_ and region III by r < r_. One copy of Taub-NUT space corresponds to the regions I, II, and III glued together diagonally, 

compare Fig. [9] 



VI. THE OBSERVABLES 



As in any other space-time orbits in Taub-NUT space-time possess some invariantly denned observables like 
perihelion shift, light deflection, the deflecton angle for flyby orbits, or the Lense-Thirring effect. These quantities 
are directly or indirectly related to the periods of the "d and the f motion. 

The two fundamental periods of the p-function where 2u>\ € K and 2w2 S C are 

c^f^U, (52) 



For a bound orbit in the region (4)_ the f-motion of the test particle is periodic f (7 + Wf) = ^(7), in particular it 
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III III 




FIG. 13: The CEO (red) and the CBO (blue) in the Kruskal-like analytic extension of Taub-NUT space-time. In the Misner-Taub 

interpretation of Taub-NUT space— time these orbits arc incomplete at the horizons. In the Bonnor— Manko— Ruiz interpretation of 
Taub-NUT space— time the orbits are incomplete at the singular axis. The orbits then consist of incomplete pieces: 3 pieces for the CEO 
(red): si — > SB, 8& — > V, — ► @, which are incomplete at SS and < g; an infinite number of pieces for the CBO (blue): 80 — > SS, SS — > , 
c € — > <2) , etc., which arc incomplete at s$ ' , SS, % ', 2> , etc. Also shown is the orbit Do (green). 



oscillates between r m in and r max with a period w? = 2uii 



UJp 



R 



dy 



2wi , 



where e% and ei are the zeros of Pziy) related to r m j n and r max . The corresponding orbital frequency is 
The ^-period of a bound orbit in the region (4)_ is simply given by 

( ) = ^ ~ ' 



(53) 



lj$ = 2 — = -2 

,/i9 min V6 J£ 



2vr 



'6 f 



(54) 
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and the corresponding frequency by £ ?L . 

The secular rates at which the angle ip and the time t accumulate are given by (for L 1 > 2Eh): 



2 



(I- 



2tt 



-a . 



(55) 



and 



r = A 

UJf 



£_t^L + Jl / (L'-2n£costf) 



A r 



Ah 2 E 



R 

2fo 



2n(cosi? + C) dd 
sin 2 1? 



(56) 



where If defined in Eq. (j49|) is evaluated at 7 Gi corresponding to the root ej, i = 1,2. The orbital frequences f2 r , £2^ 
and Op are then given by: 



Wj= 1 



2tt1 
"(9 — — p > 



V 

iip — p • 



(57) 



As discussed in [2J, \2l 
Lense-Thirring effect 



the differences between these orbital frequences are related to the perihelion shift and the 

— 2tt* 



^perihelion 



^Lcnse— Thirring 







(58) 
(59) 



Interestingly, there is no Lense-Thirring effect in the Taub-NUT space-times. The geometrical reason for this is 
the property that the orbits lie on cones. These cones are fixed in space, and so the orbit reaches the maximal value 
of d after each full revolution of <p by 2-7T. Consequently, the frequencies il v and f2,j coincide. Thus in Taub-NUT 
space-times the fixed cone plays the role of the fixed plane of Schwarzschild space-times. We expect that in the case 
of a Kerr- Taub-NUT space-time the orbital cone in the far field region precesses like the orbital plane precesses in 
Kerr space-time. 

Further possible observables in the Taub-NUT space-time are given by the deflection angle for escape and transit 
orbits for massive test particles and for light. 



VII. CONCLUSIONS AND OUTLOOK 



In this paper we presented the analytic solution of the geodesic equation in Taub-NUT space-times in terms of 
the Weierstrafi p, a and £ functions. The derived orbits depend on the particle's energy, angular momentum, Carter 
constant and on the parameters of the gravitating source. We discussed the general structure of the orbits and gave 
a complete classification of their types. 

We also addressed the properties of the geodesies in the light of the interpretations of the Taub-NUT metric following 
on the one hand Misner Q and on the other hand Bonnor Q. According to Misner and Taub [l[ one encounters 
incomplete geodesies, when the affine parameter suddenly terminates at the attempt to cross a horizon a second 
time. In their approach Misner and Taub used Eddington-Finkelstein transformations. Such kind of transformations 
lead to incomplete geodesies also in the Reissner-Norstrom space-time. This is due to the shortcoming of these 
coordinates to not provide the complete analytical extension of the space-time. In such a Kruskal-like extension of 
the Taub-NUT space-time as performed by Miller, Kruskal and Godfrey @ the incomplete geodesic behaviour at 
the horizons is eliminated. However, other geodesies are then incomplete in the analytically extended Taub-NUT 
space-time, because the singularity at the •& = 0,7r axis is retained. A periodic identification of the time coordinate 
would violate the Hausdorff property of the manifold and is thus prohibited. Following Bonnor [8j and Manko and 
Ruiz Q the singularity on the axis must then be considered as a physical singularity. 
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(a) Geodesic incompleteness of a single copy of RN 
space-time for ip+ 



(b) Geodesic incompleteness of a single copy of RN 
space— time for i/>_ 



FIG. 14: Reissner-Nordstrom space— time: ip+ and for rj = 0.4, L — 2 and i?- — 0.95. 



Appendix A: Incomplete behavior of geodesies in the Reissner-Nordstrom space— times 

The f- and ^-equations in the Reissner-Nordstrom space-time in four dimensions take the form (cf. Eqs. (O, (jTUJ) ) 



-7- = Rrn ■ 

dt _ Ef A 
g?7 Ar N 



Rm 



^f 4 -A RN (^+^) , 



A RN = r 2 



f + Tj 



(Al) 
(A2) 



where 77 is the normalized charge, 77 = and the Mino (2pj time 7 is used where f 2 c?7 = df . Solutions of these equa- 
tions can be derived by the same scheme as employed in Eq. with solution (|36|) and in Eq. ([TO)) with solution (|44| . 
(Note that in the Reissner-Nordstrom space-time the motion can be restricted to the equatorial plane and only the 
if (7) part survives.) 

The incomplete behaviour of geodesies in the Reissner-Nordstrom space-time is similar to that in the Taub-NUT 
space-time. Although the orbits are completely regular from the f— , and ^-equations, this interesting feature 



reveals itself when one investigates the ^-equation. The Eddington-Finkelstein transformations ip± — t ± 



-df 



^RN 



convert the metric into a nonsingular form. The affine parameter 7 stops to increase at some point while tp diverges 
there. This happens at the attempt to cross one of the horizons a second time, as illustrated in Fig. 1141 The similarity 
with Fig. ITTI for the Taub-NUT space-time is immediately recognized. 

The Kruskal-like extension of the Reissner-Nordstrom space-time is better suited for the description of the 
geodesies. These coordinates present an analytical extension of the space-time to the entire range of possible radial 
and time coordinates, and yield an infinite set of copies of the original Reissner-Nordstrom space-time. Consequently, 
the geodesies which seemed to be incomplete at the horizons in Eddington-Finkelstein coordinates can be smoothly 
continued through the horizons to the glued regions I, II, III of the copies. (The Carter-Penrose diagram of the 
Kruskal extension for the Reissner-Nordstrom space-time can be found, for instance, in [2H [27j.) 



Appendix B: Integration of elliptic integrals of the third kind 



We consider an integral of the type ii 



dv 



This integral is of the third kind because the function 



lv in P( v ) - Pi 

fi(v) — (p(v) — pi) 1 has two simple poles v\ and V2 in a fundamental parallelogram with vertices 0, 2loi, 2lj\ + 2u>2, 
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2w2, where 2cji and 2ll>2 are fundamental periods of p(v) and p'(v). 
Consider the Laurent series for the function fi around Vi 

fi(v) = a-x,i(v — Vi)^ 1 + holomorphic part , (Bl) 

and the Taylor series of f^ 1 about Vi 

fr l (v) = p'{v i ){v-v i ) + 0(iJ t ) . (B2) 
Comparing the coefficients in the equality 1 = fi(v)fi (v) where 

1 = h{v) (p'(vi)(v - vt) + 0(v 2 )) = a- lti p'(vi) + 0(v 2 ) (B3) 

yields a_i » = — ; — r. Thus, the function fi (v) has a residue ,} ; in u,-. 

The Weierstrafi £(i>) function is an elliptic function with a simple pole in and residue 1. Then the function 
Ai = fi(v) — 2^ = i ^?Tj^y 1S an elliptic function without poles and therefore a constant (23j . which can be determined 
from /i(0) = 0. Thus, 

^ = g p^ ' (B4) 

here p'(i>2) = p'i^uij — fi) = — p'(ui). Applying now the definition of the Weierstrafi c-function J"J ((v)dv = 
loger(v) — \oga(vi n ) upon the integral Ii we get the solution 

h = f h(v)dv = -ttU ( C(«0(« - «ln) + log ) . (B5) 



Appendix C: Integration of elliptic integrals of the type I2 = f v , r^zjz^ P 

We consider the Laurent series of $2{v) and the Taylor series of / 2 ~ 1 (w) around Uj for /2(u) = o : 

/ 2 (w) = a-2,i{v - v,A~ 2 + a-i,i(w - Vi)^ 1 + holomorphic part , (CI) 
= (p'(vi)(v - Vi) + \p"{vi){v - v t f + 0(v 3 )^j 

= (p'(vi)) 2 (v - Vl f + p'{vi)p'\vi){v - v t f + 0{v A ) . (C2) 

The function /2(f) has poles of second order in V\ and V2 such that fi{vx) = p% = ,f2{v2)- Comparison of the 
coefficients in 

1 = h(v) ((p'K)) 2 (v - Vi) 2 + f/(«i)p>i)(« " v^ 3 + 0(v 4 )) 
= a- 2 ,i (P>,)) 2 + (v- Vi) \a-i ti (p'(^)) 2 + a_ 2 ,ip'(«i)p"(«i)l + 0(v 2 ) (C3) 



yields 



1 P"(vi) (rA , 

a_2 i = T 1 CL — 1 i = T • (v^4) 

' (p'K)) 2 * (p>*)) 3 



The function p(v) possesses a pole of second order in v — with residuum and the Laurent series of p begins 
with v~ 2 . Then the Laurent series of a-2,ip(v — around Vi begins with a_2.i(w — Vi)~ 2 which is similar to 
the first term in (|C1|) . The Laurent series of a_i — Uj) around w,; begins with a-i t i(v — Thus, the 

function A2 — f2(v) — ( — — — ^ v ^^ v v ^ ) has no poles and is constant [23.1 and can be calculated 

S V (p'(».)) 2 (p'(».)) 3 J 

tarn AM -0: *--£(^ + !Oai£&!iY 

i = l \(P (P(^)) / 
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Using of p(v)dv = — ((v) + C( u in) and the definition of the cr-function the integral I2 takes the form: 

2 



f-V 

h= f 2 (v)dv = A 2 (v-v in )-^2 

Jv iD i=1 



C(v - Vi ) - ((v in - v t ) + log ■ 



'(vi) a(v in - 



1 



(C5) 
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